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ANALYSIS OF TURBULENT FLOW AND HEAT TRANSFER ON A FLAT PLATE 
AT HIGH MACH NUMBERS WITH VARIABLE FLUID PROPERTIES 1 

Uy It. (i. Dkiski.f. u and A. L. JjOkffler, Jr. 


SUMMARY 

A previous analysis of turbulent heat transfer and 
flow with variable fluid properties in sm ooth passages 
is extended to flow over a flat plate at high Mach 
numbers. Velocity and temperature distributions 
are calculated for a boundary layer in which the 
effects of both frictional heating and external heat 
transfer are appreciable. The viscosity and thermal 
conductivity are assumed to vary as a power of the 
temperature, while the Prandtl number and specific 
heat are taken as constant. Skin-friction and heat- 
transfer coefficients are calculated and compared 
with the incompressible values. The relation be- 
tween boundary-layer thickness and distance along 
the plate is obtained for various Alack numbers. 
The analytical results are compared with ■representa- 
tive experimental data. 

INTRODUCTION 

The current importance of high-speed /light has 
caused much interest in research on compressible 
boundary layers. The skin friction in high Mach 
number /light constitutes a large* part of the total 
drag. Therefore, the accurate prediction of skin 
friction is desirable for the design of high-speed 
aircraft. Prediction of heat-transfer coefficients 
in high Mach number flow is also important, 
because frictional heating of the surface necessi- 
tates cooling to prevent structural failures. 

The prediction of laminar boundary layers 
from the basic equations of momentum, energy, 
and continuity has reached a high state of develop- 
ment. A considerable amount of analytical work 
on turbulent boundary layers has also been carried 
out. In the turbulent case, however, the results 


of the various analyses disagree markedly becatise 
of the different assumptions made by the various 
authors. These analyses are reviewed in refer- 
ences 1 to 3. The introduction of assumptions 
into the treatment of turbulent boundary layers 
is at present unavoidable, since solving the prob- 
lem from the instantaneous equations of momen- 
tum, energy, and continuity alone is not yet 
possible. In some respects, however, the model 
used for solving the problem might be improved. 
In nearly all the analyses, the flow is divided into 
a laminar region, where turbulence is supposed to 
be absent, and a fully turbulent region. The 
effect of variation of fluid properties on the laminar 
region is generally neglected. Measurements of 
turbulent velocity profiles indicate that consider- 
able turbulent shear exists within the so-called 
laminar layer (ref. 4), so that a more realistic 
model for the region close to the wall than that 
used in previous analyses is desirable. 

A somewhat improved treatment of the region 
close to the wall is given in references 4 to 6, 
where the effects of turbulence and of variable 
fluid properties in this region are considered. In 
tin* region away from the wall the von Karrnan 
similarity expression Inis been considered the most 
reasonable expression available (ref. 7). In ref- 
erence 8, fully developed turbulent flow and heat 
transfer in smooth passages for air with variable 
properties are analyzed, and tin* results agree well 
with experimental data. The analysis is extended 
to the entrance regions of passages and to high 
Prandtl numbers in references 9 and 10, where 
good agreement with experiment is again obtained. 
Since the analyses apply well to entrance regions, 


1 Supersedes NACA Technical Note 4262 by If. G. Delssler and A. L. Loeffler, j r „ 1958. 
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tin' assumptions made in t lie analyses should 
apply also to a compressible boundary layer. 
The analysis is extended to flow and heat transfer 
in a boundary layer at high Mach numbers in this 
paper. (Some preliminary results were presented 
in ref. 11.) The variation of properties due both 
to frictional heating and to external heat transfer 
is considered. The viscosity and thermal con- 
ductivity are assumed to vary as a power of the 
temperature, while the Prandtl number and 


specific beat arc taken as constant. 

SYMBOLS 

A 

constant 

a 

ratio of diflusivities, €/»/« 

B 

constant 

C 

constant 

C, 

friction coefficient, 2t w! P?>u\ 

C ‘ P 

specific heat of fluid at constant pres- 
sure 

D 

constant 

d 

exponent for viscosity variation with 
temperature, taken as 0.6S for aii 

E 

constant 

n 

enthalpy 

h 

h ea t - 1 ra nsf er eoeffi cie n t , q w j taw) 

k 

thermal conductivity 

u 4 

Mach number based on f roe-stream 
properties and velocity, Uil^'yRts 


constant, 0.109 

Pr 

<1 

Prandtl number, c P p!k 
energy transfer in ^-direction per unit 
time per unit area (as defined by eq. 
(A 14)) 

It 

j)(‘rfeet gas constant 

Re, 

llevnolds number based ou x, xuspups 

R<« 

Reynolds number based on 0, Bu^pbipb 

St 

Stanton number, hlc v Ubpb 

T 

total temperature, t J r{u 2 j2v P ), deg abs 

T + 

t o t al- temp er at ure p ara meter , 

( t T)C v T w 1 — (T/t „,) 

(jw'V T wl Pw ^ 

t 

static temperature, deg abs 

t+ 

t cm porat ure paramo ter, 

(t w —t)C p T u 1 — (t/t„) 

q^Twjpu, P 

t+' 

temperature parameter, 
2(t,r— t)c v p„ 1 — (t/tj) 

T w OC 


n. 

r 

x 

V 


yt 


P 

7 

5 

5+ 


5/1 

5/1 


€ 

V 

e 


o + 


velocity in ^-direction 
velocity parameter, u/\ tJp w 
velocity in ^/-direction 
longitudinal distance along plate 
distance perpendicular from plate 

„ .. j y \ r w!ptp 

v ail distance parametei, 

Pwi P w 

lowest value of ?/ + for which equation 

for region away from wall applies 

frictional-heating parameter, T w j2c p t w p w 

. c , Qw\ t wIPid 

heat-transfer parameter, — — 

1 pi tr ' w 

ratio of specific heats, taken as 1.400 
for air 

flow boundary -layer thickness 

flow boundary-layer-t hickness param- 

Her, iVjs/fh? 

P wj P to 

tiiermal boundary-layer thickness 
1 1 1 er m al bound ary-1 aver- tliieki less pa- 

5/1 \ T w! Pic 
Pic I Pic 

eddy diffusivity of momentum 
eddy diffusivity of heat 


rameter, 


t air 

V 


t e mpera ture -reeo ve rv factor, 
momentum thickness, 

r p 

t /n Pi ^ 5 V ^'6/ 

mo m en t um- 1 1 1 i ck i ) ess para me ter, 


l\ZtJ_Pjc 
P id Put 

constant, 0.30 
viscosity 
density 

shear stress, force per unit area 

Subscripts: 

aw pertaining to adiabatic wall conditions 

i incompressible; constant fluid proper- 

ties 

pertaining to wall 

pertaining to edge of boundary layer 
or free stream 

pertaining to edge of wall layer 
Superscripts: 

* reference 

' pertaining to fluctuations from time 

average except in t +/ 
time-averaged value 


w 

5 
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ANALYSIS AND DISCUSSION 


BASIC EQUATIONS 

The instantaneous velocities, temperature, and 
fluid properties in the equations of momentum, 
energy, and continuity can be divided into mean 
and fluctuating components. If time averages 
are taken, the following equations for shear 
stress and energy transfer, applicable to flow in 
a boundary layer, are obtained (appendix A): 


q= — k 


(\u —7 

T=jLi “j p U 

(1 V 

1/ 


a) 


< 1 ?/ 



pc p t'v' — up 

d y 

+ up u'v' 

(?) 

specific heat 

is 

assumed. 

The 


iU , — 

i\y 


bars denote time averages, and the primes indicate 
fluctuating components. Equations (1) and (2) 
are the same as equations (A 9 ) and (Al4) in 
appendix A if the bars over the time-averaged 
velocities, temperatures, and properties are 
dropped. The various terms in equations (l) 
and (2) may be interpreted as follows: 


Au 



molecular shear stress 


— P u'v' turbulent shear stress 


—k molecular heat transfer 

dy 

pc p t'v' turbulent heat transfer 


— vu 4^- molecular shear 

<* y 


work 


up u'v' turbulent shear work 


r=(/i+pt)^ ( 3 ) 

( I~ — (k+pCpth) u(p-\-pt) ^ ( 4 ) 

The physical significance of e and e„ lies in the 
fact thate/(p/p) is the ratio of turbulent to molec- 
ular shear stress (ref. 12), and e h /(k/pc p ) is the 
ratio of turbulent to molecular heat transfer. 
Equations (1$) and ( 4 ) can be written in dimension- 
less form as 


P € \ d ?/ + 

r tr Pw PwjPw' d?/ + 

and 


JL = (JL 1 j_ JL a € \ + 

(]w \k IQ I t W Pw P-U'ipw) d?/ + 

_ 2 - u+ ( M + p - — *_Y— - 
P \Pw Pw PwjPw/^y^ 


( 6 ) 


The subscripts w refer to values at ?/=(); that is, 
at the wall. The quantity a is a frictional- 
heating parameter that is an indication of the 
variation of properties due to frictional heating, 
and p is a heat-flux parameter that is an indication 
of the variation of properties due to heat transfer. 
The parameter a is always positive or zero, a 
value of zero characterizing low-speed flow (i.e., 
A/j= 0). A zero value of ft refers to a vanishingly 
small heat transfer or an insulated plate. A 
positive value of p indicates heating of the fluid, 
while negative p means that the fluid is being 
cooled. It is sometimes convenient to write 
equation (6) in the following alternative dimen- 
sionless form : 


Equations (1) and (2) suggest the form of the 
turbulent transfer equations but contain the un- 
known quantities u'v' and t'v', so that assumptions 
must be made before solutions can be obtained. 
For making these assumptions it is convenient 
to introduce the relations 

— — dw , -7—7 lit 

u v ~ — and t v - 

(l?/ d y 

where e and are the eddy diffusivities for 
momentum and heat transfer, the values of which 
depend upon the amount and kind of turbulent 
mixing at a point. When these relations are 
introduced, equations (1) and (2) become 


P q/ k 1 p ^ j_\ dt [' 

^ q w \k ip 1 i ip p w p, wj P w d// + 

- 2 u+ ym ( 7 ) 

\pir pw Pw/Pw/ d y 

This equation is particularly convenient when 
>3 = 0 , for which case equation ((>) becomes indeter- 
minate. 

EXPRESSIONS FOR EI)I>Y DIFFUSIVITY 

In order to make practical use of equations (5) 
to ( 7 ), the eddy diffusivity e must be evaluated 
for each portion of the flow. For this purpose the 
boundary layer is divided into two portions termed 
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the “region away from the wall” and the “region 
close to the wall.” 

Region away from wall. — In the region away 
from the wall, it is assumed that the turbulence 
at a point is a function mainly of local conditions — 
that is, of the relative velocities in the vicinity of 
the point (ref. 13). This is probably not a good 
assumption in the region near the edge of the 
boundary layer, where considerable diffusion of 
the turbulence occurs (ref. 14) and, in addition, 
the turbulence is intermittent. However, in that 
outer region the velocity or temperature gradients 
are so small with respect to these gradients nearer 
the wall that the error in calculated velocities or 
temperatures should not be large. A Taylor series 
expansion for u as a function of transverse dis- 
tance, then, indicates that t is a function of 
du/d?/, d 2 u/d?/ 2 , d 3 u/d?/ 3 , and so forth. If, as a 
first approximation, e is considered as a function 
onlv of t lie first and second derivatives, and dimen- 
sional analysis is applied, 



This expression was obtained by von K firman and 
is generally known as the Ivarman similarity 
hypothesis (ref. 7). The constant k is to be deter- 
mined experimentally. 

Region close to wall. — In the region close to the 
wall it is assumed that € is a function only of 
quantities measured relative to the wall that is 
of u and yr This assumption includes, to a first 
approximation, an effect of the derivative du/i\y. 
Since the flow becomes very nearly laminar as the 
wall is approached, the first derivative approaches 
the value ujy and hence may be omitted, since u 
and y already appear in the functional relation. 
By using dimensional analysis, 

e~e(u,y)—n 2 uy (9) 

whore n is an experimental constant. 

Equations (8) and (9) can be considered as 
reasonable tirst approximations for e. Whether 
these approximations are adequate or not can at 
present he determined only by experiment. 

- Kofert'tuv 10 sli«ws that the kinematic viscosity has an effect on « in the 
region very float! to the wall. However, that effect becomes important only 
for heat or mass transfer at I*ran<i tl or Schmidt numbers appreciably greater 
than l. 


Determination of experimental constants. — The 

constants n and k were determined from pipe data 
in which the properties were essentially constant. 
Equation (5), with equation (8) or (9), was inte- 
grated (constant properties and r) for the regions 
close to and away from the wall in reference 4. 
The molecular shear stress was neglected in the 
region away from the wall, and the well-known 
K arman-Prandtl logarithmic equation was ob- 
tained in that region. In matching the two 
solutions it was assumed that the velocity is con- 
tinuous at the junction of the two regions. 

The integrated equations (ref. 4) for the regions 
close to and away from the wall are plotted in 
figure 1 with the constants n — 0.109 and k= 0.36 
determined from pipe data (refs. 4 and 14). The 
data indicate that the equation for the region close 
to the wall applies for ?/ + <^26, and the equation 
for the region away from the wall applies for 
2/+^>2f>. Included in the plot are data for a low- 
speed boundary layer with zero pressure gradient 
from reference 15. The agreement with the curve 
is satisfactory. 

The values for the constants n~ 0.109 and 
k =: 0.3G should apply to flow with variable as well 
as constant properties if the basic assumptions 
mad k for e in the preceding sections apply to 
variable properties; that is, if € — e{u,y) close to 
the vail and e = e(d?//d//, (Vujdy 2 ) away from the 
wall The constant t v|, however, requires fur- 
ther consideration and is discussed in the next 
section. 

ADDITIONAL ASSUMPTIONS 

Ii addition to the assumptions for eddy difFti- 
sivity discussed in the preceding section, several 
additional assumptions must be made for solving 
equations (5) to (7). 

Variation of properties with temperature. — For 

gasts, the viscosity varies approximately as t d , 
where d has an average value of 0.08 for tempera- 
tures between 0° and 2000° F. The Prandtl 
number (Pr= 0.73) and specific heat c p arc as- 
sumed constant, because their variations with 
temperature are of a lower order of magnit ude t han 
the variations of the other properties. If c p and 
Pra id 1 1 number are considered constant, the 
thermal conductivity A* will vary with temperature 
in the same way as the viscosity, or as t d . For 
constant pressure across the boundary layer, the 
density p is inversely proportional to t. 
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Figure L — Generalized computed velocity distribution for constant-property turbulent flow (ref. 4) with constants in 
equations determined from pipe data compared with data for low-speed boundary-layer flow on a flat plate. 


With the preceding assumptions, the property 
ratios in equations (5) to (7) can be written as 


Pw & w \tw/ \tw/ 


P _ 1 

Pw ^/itc 

From the definitions of (3 and t+, 
t 


i t to 


= 1 -/W+ 


( 10 ) 

( 11 ) 

( 12 ) 


or, if equation (7) rather tlmn equation (6) is used, 
t 


t 


-=1 -at + f 


(13) 


The property ratios in equations (5) to (7) can 
therefore be written in terms of (3 and t + or a and 
t + '. 

Variations of r and q across boundary layer. — 

The momentum equation (A7) indicates that, for 
a flat plate (zero pressure gradient), 
dr/dy— d(p.du/dy)!dy=() at t lie wall. Since r is 
zero at the edge of the boundary layer, the actual 
variation of r across the houndary layer might be 


expected, in general, to lie between a linear varia- 
tion (r/r M — 1 — (y/5)) and t/t w =1. Data on low- 
speed isothermal flow over a flat plate (ref. 15) 
show that this type of variation does exist, except 
in a narrow region near the edge of the boundary 
layer. For determining the sensitivity of the 
velocity or temperature profile to shear-stress 
variation, it should therefore be sufficient to com- 
pare the profiles for a constant and for a linearly 
varying shear stress. Appendix B shows that 
rlT w —qlq w for a flat plate if the Prandtl number is 1. 

Figure 2 shows u + or T + plotted against y + for 
a Prandtl number of 1 for both a constant and a 
linearly varying shear stress and energy transfer, 
where T + is the total-temperature parameter. 
Curves are shown for (3=0 and a=0, 0.003, and 
0.008, which cover much of the range of Mach 
number and Reynolds number of interest. The 
equations for calculating the curves are given in 
appendix D. The equation for the region away 
from the wall was taken to apply for ?/ + >30 
rather than >26 when the shear stress was 
variable, in order to give better agreement with the 
data for constant properties. The curve for 
a — 0.008 is cut off at the point shown because the 
Mach number becomes infinite, as can be seen 



Velocity parameter, u+, or total - temperoture parameter 
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from equation (D3) (for otu +2 = 1). The curves 
indicate that variable shear stress and energy trans- 
fer have but a slight effect on the velocity and 
temperature profiles. Similar curves were ob- 
tained in figure 1 1 of reference 6 for 09 *i) and or— 0. 
The same conclusions should apply to Fraud tl 
numbers differing slightly from 1, so that the 
effects of the variations of r and q across the 
boundary layer are neglected for solving equations 
(5) to (7). 

Ranges of applicability of equations for flow 
close to and away from wall. — It was determined 
from the data for constant properties that the 
lowest value of ?/ + for which the equation for the 
region away from the wall applies is ~2 6 when 

the variation of shear stress with y is neglected 
and the molecular shear stress is neglected in the 
region away from the wall. The question arises 
as to how yt varies when the properties are vari- 
able. The simplest assumption is that yt is 
constant and equal to 26. This assumption, which 
implies that the wall properties govern the thick- 
ness of the wall layer (yt is 

similar to von Karman’s assumption (ref. 16). 
Figure 12 of reference 6 shows that essentially the 
same curves are obtained when the molecular 
shear stress is neglected in the region away from 
the wall as when it is considered, the difference 
being that, when the molecular shear stress is 
included, yt has the constant value, of 16 rather 
than 26. 

Another assumption, which might be somewhat 
more reasonable than assuming yt constant, is 
that yt occurs at a given constant ratio of turbu- 
lent to molecular shear stress ej(pfp). That is, 
the turbulence changes from that described bv 
equation (9) to that described by equation (8) 
when the ratio of turbulent to molecular shear 
stress reaches a certain value. In this case the 
more complete equations are used for the region 
away from the wall, in which the molecular shear- 
stress and heat-transfer terms are retained and 
the slopes of the equations for flow close to and 
away from the wall are matched at yt (^=16 
for f$=a~ 0). 

In figure 3, or T + is plotted against ?/ + for a 
Prandtl number of 1 using the two assumptions 
for yt discussed in the preceding paragraphs. 
Curves are shown for 0 and a = 0, 0.003, and 
0.008. The equations for calculating the curves 
are given in appendix C. The curves indicate 

n<;s;U:V 00 -- 2 


that the velocity and temperature profiles are 
apparently insensitive to the assumption used 
for yt . Similar results were obtained in figure 
13 of reference 6 for (3^ 0 and a— 0. The simpler 
procedure of neglecting the molecular shear stress 
and heat transfer in the region away from the 
wall and assuming yt = constant = 26 is therefore 
adopted in the following calculations. 

Ratio of eddy diffusivities for heat and momen- 
tum transfer. — In most analyses the ratio of eddy 
diffusivities a that occurs in equations (6) and 
(7) is set equal to 1; that assumption has given 
heat-transfer coefficients in good agreement with 
experiment (ref. 8). It is of interest that FrandtFs 
mixing-length theory, which assumes that a 
turbulent particle moves a given distance and 
then suddenly mixes with the fluid and transfers 
its heat and momentum, gives a value of a — 1 . 
Although the actual turbulence mechanism may 
be more complicated than indicated by that 
theory, it does indicate that a value of a on the 
order of ] is not unreasonable. 

In the present analysis the assumption of 
a — 1 is retained, but in some cases the calcula- 
tions are also carried out for a— 1.07 in order to 
determine the effect of varying a, A ratio of 
diffusivities of 1.07 was obtained from some pre- 
liminary experiments on recovery factors for 
fully developed flow in a tube. 


VELOCITY AND TEMPERATURE DISTRIBUTIONS IN 
BOUNDARY LAYERS 

For obtaining velocity and temperature dis- 
tributions close to the wall, equations (9) to (13) 
are substituted into equations (5) to (7). Equa- 
tions (5) and (6) become, in integral form, with 

t/t u> y_! — l } 


r y + 

Jo 


d?/ 4 


(1 -Pt + ) d - 


t + = 


-r. «= 


i 

(l+2"u+)<b/ 


+ n 2 u + y + 


Pt + ) , fl 

/V 1 1 07+™^ 


(14) 


(15) 


.ii+ 


Equations (14) and (15) can be solved simul- 
taneously by iteration; that is, assumed relations 
between u + and ?/ + and and ?/ + are substituted 
into the right sides of the equations, and new 
values of u ¥ and are calculated by numerical 
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integration. These new values are then sub- 
stituted into the right sides of the equations and 
the process is repeated until the values of u + and 
do not change appreciably. Equations (14) 
and (15) give the relations among u + , t+, and y + 
for various values of a and 0 for flow close to the 
wall (?/ + <2(5). For 0—0 and 0, t + becomes 
infinite, so that equations (5) and (7) must be 
used. These equations, with equation (111), 
become 




nhi + y + 


(16) 



g+2^)dy‘ 


(1 -at +r ) d . 1 2+ + 

Pr -+T^i^ nUy 


(17) 


Equations (16) and (17) are solved similarly to 
equations (14) and (15). 

In the region away from the wall, the molecular 


shear stress and molecular heat transfer are 
neglected. Dividing equation (6) by equation 
(5) gives, with t/t» — g/g„= 1, 

1+2 i“ + = a d £ (18) 

Integrating equation (18) from yt to y + gives 

t+ = tt+—--t + £. u +2 “ ut* (19) 

a a 0a 0a 

From equations (11), (12), and (19), 


i-ptt+^ u ' +- ur 2 — u+ -- u+ 2 
a a a a 


Substitution of equations (8) and (20) into (5) and 
one integration give, for the region away from 
the wall, 

k _i I " 2 ortt ++/3 — J 

jr <!>/,+ Lv^ 2 + 4a < a^a6<t+l> u t+ au l 2 ) J A) , \ 


K ^r e 



(a) Velocity distribution. 

Una kk 1. Predicted generalized velocity and temperature distributions for air with Heat transfer and frictional heating. 

Prandtl number, 0.71. 
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By lotting 

2 =—= sin -1 r - 

V“/» Lv 

and integrating equation (21), 


2tm+-f/S _ 
0 2 +4a(n — afitt +put +«m 


— ( 22 ) 

i +2 )J 


r 


AV 


1 + 


( \a/a \ a/a . \a/a \ , T , _ s 

(™ S A" 2+ « (2») 


a k z 


whore 2 is given by equation (22) and z l is the 
value of 0 at yt = 20. Equations (22) and (24) 
give the relation between u + and y + for various 
values of a and 0. The quantity /+ can then be 
calculated from equation (19). 

Por = 0 and 0, l)ecomes infinite, and 
t+ f =(0/a)t + must be used. Equation (19) be- 
comes, in terms of t + 


1 he constant K is evaluated in the usual way 
by letting dw + /d ?/ + — 00 at y+=0 in equation (21) 
(ref. 7) and substituting (21) into (22) and (23) 
at y+ = 0. 3 Bv using this procedure, K= 0. To 
determine K t) , set u+=ut when y + =yf, Then, 



3 This assumption can he avoided by including the molecular shear stress 
and heat transfer in the region away from the wall and evaluating K by 
assuming a continuous velocity derivative at yf (fig. 12, ref. (>). This assump- 
tion gives essentially t lie same results as that made iu the text. 


t +, =ty 


£ 

aa 


u\ 


uV 0 

a T aa 


u + + 


ir 


( 25 ) 


Equations (22) and (24) apply to the case for 
0=0 and a 9^0 if atf is substituted for 0tt in 
equation (22). 

For a=0, equation (23) becomes indeterminate, 
and equation (21) for zero frictional heating from 
reference 6 can be used. 

Typical velocity and temperature distributions 
for various values of the frictional-heating param- 
eter a and of the heat-flux parameter 0 are 
presented in figures 4 and 5. Positive values of 
0 correspond to heat addition to the air; negative 



(b) Total-temperature distribution. 

Fjoriiu 4.— Concluded. Predicted generalized velocity and temperature distributions for air with heat transfer and 

frictional heating. Prandtl number, 0.73. 
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values, to heat extraction. The curves of u + 
against ?/ f (tig. 4(a)) indicate considerable flat- 
tening of the velocity profile as either a or 0 
increases positively. This is caused by the 
decreasing temperatures in the outer regions of 
the boundary layer compared with the wall 
temperature when either the Much number is 
high (high a) or the heat transfer from the surface 
to the air is high. Thus, the density is higher 
in the outer regions of the boundary layer, w itli 
consequent flattening of the profile (eq. (i>)). 
Negative values of (3 produce the opposite effect. 
For certain combinations of a and 0 (with 0 
negative), the effect of 0 on the curves should 
tend to cancel the effect of a, and the resulting 
profile should not differ greatly from the a— 0=0 
curve. The curve for a= 0.002, 0— — 0.05 in 
figure 4(a) is close to the curve for <x=0=O. 
Included in figure 4(a) for comparison are experi- 
mental data from reference 17 for an a of 0.00176, 
0=0, and a corresponding Mach number Ms of 
2.82. The data arc in reasonable agreement with 
the predicted profiles. 


In figure 4(b), is plotted against // + for 
various values of a and 0. The total -tempera lure 
parameter T + is plotted rather than f*”, because 
the trends are somewhat more consistent, although 
some crossing over of the curves occurs even 
with 7 ,+ . For calculation purposes, a better 
representation can be obtained by plotting 7 + 
againrt u + . The quantity 7 M is related to t' by 
the relation 

T + = t + —%u +t (26) 

a 


SKIN -FRICTION COEFFICIENTS 

Tli • skin -frict ion coefficient is defined as 
/< _ 

< / o 

PbU\ 


(27) 


where the subscript 8 refers to values outside the 
boundary layer. Equation (27) becomes, in 
dimensionless form, 

2a-atr_)jM[-0fr) (28) 

Us 2 Us 2 

For comparison with experimental data, it is 


fi _ 

f Ps ut 2 ~ 
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convenient to introduce the momentum thickness, 

l--)dy (29) 

Jo psUi\ Us/ J 

which in dimensionless form is 

r-d-.inj' <:,0) 

Then the Reynolds number based on the momen- 
tum thickness and free-stream properties is 

Re, = 6,1 iPi =e+ni - * ~ (3 1 ) 

M<5 M<5 Ptc 

where the property ratios are obtained from 
equations (10) to (13). The Mach number for a 
perfect gas is 


which in dimensionless form is 


V w-Wh/o ^ 

If values are given to a, 0, and 5 + , where is 


the value of at the edge of the flow boundary 
layer, then values of ?/ 8 + , / g + , and so forth can be 
read from curves similar to those in figures 4 and 
5. Values of C f} lie$, and Ms can then be calcu- 
lated from equations (28), (31), and (32). This 
procedure assumes that the thermal and flow 
boundary layers are of equal thickness. From 
the calculations in a later section, where relations 
between boundary-layer thickness and distance 
along the plate arc calculated, it can be shown 
that this is a good assumption for gases when the 
thermal and flow boundary layers begin at the 
same point. For the case of Pr~a— 1 the assump- 
tion holds exactly, as can he seen by substituting 
?/ + — 7 T+ into equations (42) and (43), which are 
then identical. 

Predicted skin-friction coefficients are plotted 
against life in figure 0(a) for various values of 



(a) Insulated plate; heat-flux parameter fi, 0. 

Figure 6 . — Variation of predieted skin-friction coefficient with momentum-thickness Reynolds number and Mach number 

and comparison with experiment. Prandtl number, 0.73. 
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(b) i„ o r*. 

Kici'ui-: (>. ( ’onrludtMl. Variation of pmliettd skin-frict ion (.‘oWJicion with moment um-thieknvss Reynolds number and 

Mm*1i number and comparison with experiment. Pram.lt l number. 0.73. 


Mach number for an insulated plate (fi 0). 
'These curves are for €/,/«-- a 1 . Tlie effect on the 
curves of changing a to 1.07 was negligible. The 
values of (' f decrease considerably as Mach 
number increases. Included in the plot are 
experimental data of a number of investigators 
for Mach numbers up to 4.93. In general, the 
data are in good agreement with the predicted 
curves. 

The ratio of the friction coefficient- to the in- 
compressible coefficient is plotted against Mach 
number for various values of Iiee for fi— Q in figure 
7(a). The values of C f \ (' fi decrease with AY*, but 
at a decreasing rate. For comparison purposes 
the analytical curve of ('fj(\,i against M 0 for a 
value of life of <>()()() is plotted in figure 7(b) 
together with data taken near this value of lies. 

If heat transfer occurs between the plate and 
the stream, it is convenient to specify the ratio of 
the actual wall temperature to the adiabatic wall 
temperature for a given Mach number and Ile$. 


For an insulated [date the adiabatic wall temper- 
ature r lav be written as 

* (33a) 

2(i P 

where rj is the temperature-recovery factor, the 
calculation of which is discussed in the next 
section, Equation (33a) can be- written in dimen- 
sionless form as 

~--=. ' (:«»» 

faic 1 ““ +VOiU d ‘ 

Figure 0(b) is similar to figure 0(a), except that 

the ph te is now cooled (b r /f-«io = =b.* f 5)- The trends 
are sii lilar to those of figure 0(a), but all the 
curves are displaced upward. This increase in 
friction coefficient was also obtained for flow in a 
tube with cooling (ref. 5). Also included in this 
figure ire wind-tunnel data (ref. 18) obtained at 
high Mach number using nitrogen as the working 
fluid. The agreement with theory appears to be 
within experimental error. 
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with properties evaluated at the frcc-stream tem- 
perature, is given by 



CpUbPs . + + P 5 

h u s 

Pw 



1 -fl-icMt* 


(35) 


where equation (33a) is used. The temperature 
and density ratios are determined from equations 
(11) and (12) or (13). For 0=0 and equation 
(35) becomes indeterminate. For that case set 
1 — (Ult 1t )== atf' and tj — (a/0)#'. With these 
substitutions equation (35) becomes 

Sl= —r 7 . A, +2 ; (36) 

aW a(p»/p«-)(*« —nuV) 

From equations (7) and (5) (q/q w — T /T tr =l) f 


r~- 

4 a I k J 

1/ 0 ft ID ^ 1 


dr 


1 P 6 

H — — f — 

1 w Pw M tr/ P m; 


+ f\ 07) 

V 0 » tc i l IP P w P If/ P If 


The second integral in this equation can be re- 
placed by (#'), 9 = 0 . But equation (33a) can be 
written in dimensionless form as 


yfot* — ( t 1 )(3=o (38) 

Substituting equations (37) and (38) into (36), 
with the second integral in equation (37) replaced 
by (#')*-o, gives 



1 

<iff+ 

tip € 

/. /j„ + a 

if 1 f w P ic Pwj P tr 


(39) 


For evaluating equation (39) in the region close 
to the wall, tj{p lc jp tc ) — n 2 u + y + . For the region 
away from the wall, e could be obtained from 
equation (8). However, it is more convenient to 
obtain e from equation (5), which for the region 
away from the wall becomes 


P_ y __ 1 _ 

PwPw/Pw du+fi\y+ 


culated from equations (39), (32), (33b), (30), (31), 
and (10) to (13). 

Predicted Stanton numbers are plotted against 
Rfe for various Mach numbers for tjt a w~ 1 in 
figure 8(a). The case of t w /t a10 = 1 is a limiting 
case that can be approached as closely as desired 
by making the heat flux small. When t w ft aw — 1, 
there is no effect of variable properties due to 
heat flux. The Stanton numbers in figure 8(a) 
show trends similar to those of the friction coeffi- 
cients in figure G(a). Included in the figure are 
experimental data for low heat flux obtained by 
a number of investigators. In general, the data 
are in good agreement with the predicted curves. 4 
The curves in figure 8(a) are for eje=a = 1. 
Similar curves for an a of 1.07 were 3 to 5 percent 
higher for a Mach number of 0, but the difference 
decreased at higher Mach numbers. The curves 
for o — l are in slightly better agreement with the 
data than those for a = 1.07. 

The ratio of Stanton number to the incompress- 
ible Stanton number is plotted against Mach num- 
ber for various values of Re e for t u ft ati} =l in figure 
9(a). These curves are very nearly the same as 
those for C f jC fii in figure 7(a). 

Figure 8(b) is similar to figure 8(a), except that 
twitaw = 0-5. As was the case for the friction 
coefficients in figure 6(b), the Stanton numbers 
increase as t w /t aJ o decreases. The corresponding 
plot of St/Sti against Mach number for various 
values of R/ j e and t w jt aw = 0.5 is shown in figure 

9(b). 

Temperature-recovery factors, as calculated 
from equation (38), are shown in figure 10 for 
Mach numbers from 0 to 8. Curves are shown for 
a— 1 and er=1.07. The curves for a— 1.07 are in 
somewhat better agreement with most of the 
experimental data than those for a= 1. This 
does not mean that an a of 1.07 should be used 
for calculating heat transfer or Stanton numbers. 
According to Reichardt’s hypothesis, the value of 
a should be close to 1 at the wall and increase 
with distance*, from t he wall (ref. 12). The 
temperature profiles for calculating heat-transfer 
coefficients are very steep near the wall, so that 
the important changes with distance take place 
near the wall whore a is close to 1. In the case 


Equation (39) can be used for 0 = 0 or 0=^0. For 
given values of a, and 0, values of Stanton 
number, Mach number, t lc /t aiCi and Re 0 can be cal- 


4 In a recent work, Triiss (ref. 34) extended the present analysis to the case 
of mass transfer for a Schmidt number of 2.5. The results were in very good 
agreement with experimental data for sublimation mass transfer of naph- 
thalene from a flat plate to air at Much number from 0.5 to 3.5. 


nnsai r»— no 
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(a) Insulated plate; heat-flux pa ameter /3, 0. 

t'icrKK 9. — Variation of St/Stj with Mach number for various valutas ol moment inn-thickncss Reynolds number. Prandtl 

number, 0.73. 


of recovery factors, however, the plate is insulated, 
so that the temperature gradient is zero at the 
wall. The gradients near the wall will therefore 
he smaller than in the ease of heat transfer, and 
important changes of temperature with distance 
might occur in regions away from the wall where 
a is somewhat greater than 1. 

Figures 11(a) and (b) show the curves of the 
Reynolds analogy factor 2StjC f against Mach 
number for various Re$, for iu/iatc equal to 1.0 
and 0.5, respectively. If Reynolds' analogy held 
strictly (Pr=a—- 1), the Stanton number would 
be equal to one-half the friction factor and 
2St/C f would be unity. Figure 11 shows a varia- 
tion of the Reynolds analogy factor over the range 
1.005 to 1.280. In general, 2 St/C; increases with 
increasing Mach number and with decreasing 
Reynolds number and increases slightly with 
decreasing t u Jt alr at the higher Mach numbers. 
These results are in approximate agreement with 
those of Rubesin (ref. 19), who estimated that 
2St/C f would be in the range 1.18 to 1.21 at least 
up to d/5 of 5. 


To obtain approximate values of St as a function 
of Ree and M& for t w /t a w other than 1.0 and 0.5, it 
is recommended that the results of figure 11 be 
interpolated or extrapolated to give the value of 
2 St/Cf at the desired condition. Then C f can 
be found from equations (34b) and (34c) as pre- 
viously described, and thus the value of St is 
obtain d. 

RELATION BETWEEN BOUNDARY-LAYER THICKNESS AND 
DISTANCE ALONG PLATE 

From the results given in the preceding sections, 
the stin friction or heat transfer for a given 
bound try-layer or momentum thickness can be 
calculi ted. In order to calculate the relations 
between thermal or flow boundary-layer thickness 
and distance along the plate, the well-known 
integer 1 momentum and energy equations may 
lie usi (1. These equations mav be written as 
follow? for a flat plate (zero pressure gradient): 
d f 5 

T w= <Lr J ^ u *~ d0) 

d r 5 /i 

[Cvpucr- 7’ 4 )(Ii/] 

(1./ J 0 


(41) 
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Fun :hk 9. — Concluded. Variation of <S7/.S7, with Mach number for various values of momentum-thickness Reynolds 

number. Prandtl number, 0.73. 


If 5 = 8 n ~ 0 for :/;~0 and equations (40) and (41) 
are integrated with respect to x, they become, in 
dimensionless form. 


AVv 

T 1 -, 

<r- .1 

f + - p -n+(ir t - 



J 0 P u' 


_ M<5 M 5 ,./() Pic 

' J 





(42) 


r [1 

Tp 


(i r- 7 ! 



Jo 

P tr 

LMi IjJo Pu 



— T+)i\y+ (41V) 


where the bracket for the upper limit of integra- 
tion refers to the value of the variable of integra- 
tion at that point. These equations give the 
relations between <5 + and Re x and <5* and R( x . 
The property ratios are obtained from equations 
(10) to (13)/ 

Equation (42) can be written in the more 
convenient form 



d life 

(V 


(44) 


The Reynolds number based on momentum 
thickness lies is plotted against Ih x , as found from 
equation (44), for an insulated plate in figure 12(a). 
The value of R< e decreases at a given Re x as the 
Mach number increases if the free-stream proper- 
ties remain constant. This is caused (eq. (44)) 
by the decrease of friction factor with increasing 
Mach number (fig. 0(a)). Data included in 
figure 32(a) agree reasonably well with the 
analytical curves. 

Figure 12(b) is similar to figure 12(a), except 
that t w j1 a 0.5. For given values of Re x and Mach 
number the values of Rie are generally a little 
higher for t w /f tnc of 0.5 than for t w /t (l!r of 1. This 
trend can l)e understood from examination of 
equation (44), since C f is liigh(*r for t 1c /t aw of 0.5 
(fig. 6(b)) than for f rr t aw of 1.0 (fig. 0(a)). 

Predicted skin-friction coefficients for an insu- 
lated plate are plotted against Re z in figure 13(a). 
The trends with Mach number are similar to 
those obtained when C f is plotted against Re 9 


Temper ature - recovery factor 
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(a) Insulatrd plate; heat-flux parameter 0. 

(b) f w ?t a>r 0.5. 

Iiui rk 11.- Predicted variation of Reynolds analogy factor with momentum-thickness Reynolds number and Mach 

number. Prandtl number, 0.73. 


but are less pronounced, because the boundary- 
layer thickness at a given decreases with increas- 
ing Mach number. Experimental data for low- 
speed How included in the figure are in good 
agreement with the predicted cum* for a Mach 
number of zero. Data for higher Mach numbers art* 
also in reasonable agreement with the predicted 
curves but art* somewhat, more scattered than tin* 
data in figure (>(a), where (\ is plotted against 
Rf§. This scatter is apparently caused by uncer- 
tainty as to the point at which the boundary 
layer actually starts in a supersonic 4 How. 

In figure 14(a), the theoretical curves are re- 
plotted as C against Much number for 
various Reynolds numbers based on x . The 
effect on CfjCf' t of varying Re x becomes appreciable 
at high Mach numbers. 

Stanton numbers for an insulated plate are 
plotted against Re x for a Mach number of zero 
in figure 15. Curves for higher Mach numbers 


and for t w jt aw of 0.5 involve considerably more 
calculation and were not obtained. 

higures 13(b) and 14(b) are analogous to 
figures 13(a) and 14(a), respectively, except that 
they arc tor t w !t aw of 0.5. The friction factors, as 
expected, are higher for the larger rates of cooling. 

REFERENCE TEMPERATURE FOR He. RESULTS 

The customary use of a reference temperature 
concept requires that the Reynolds number 
dependence of O^and of St be tin* same for all 
Mach numbers so that ( 1 r /( V.i and Sf/Sft should 
not lx* functions of Reynolds number. Examina- 
tion of tin* predicted curves of ( -/T7.i and St/St t 
against Mach number as shown in figures 7(a), 
7(c), 0, and 14 shows, however, that ('f!( 1 f, i and 
StjSti are strong functions of Reynolds number at 
the higher Mach numbers. Therefore, the present 
theory cannot lx* represented accurately by one 
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(b) tjt aw = 0.5. 

Figure 12. — Concluded. Predicted variation of moment um-thickness Reynolds number with longitudinal-distanci 
Reynolds number and Mach number and comparison with experiment. Prandtl number, 0,73. 



Skin- friction coefficient. 
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(a) Insulated plate; heat-flux parameter /3, 0. 

(»>) / 0 . 0 . 

I'KiruK 13. Predicted variation of skin-friction coefficient with loi gitudinal-distance Reynolds number and Mach 
number and comparison with experiment. Prandtl number, 0.73. 
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(a) Insulat-rd plate; heat-llux parameter /3, 0. 

Figure M. — Variation of CfjC / ,i with Mach number for various values of longitudinal-distance Reynolds number. Prandtl 

number, 0.73. 


reference temperature valid for all Reynolds 
numbers. 

For purposes of comparison, however, l lie 
results obtained by using Eckert T s suggested 
reference temperature (ref. 3) are shown in 
figure 14. Agreement with the present theory 
for Re x — 10 s is quite good. Use of Eckert's 
reference temperature method to solve for ('/!( V,t 
is recommended, then, if large values of Re x (near 
10 s ) are considered. In order to solve for ( /, t he 
value of (Jf t i can be taken from the curve for 
M& = 0 in figure 13(a). An approximation (within 
5 percent) to 1 his case is 

6/— 0.0292 /!e x - o m 

In order to solve for values of C f at values of 
tjtaw other than 1.0 and 0.5 and for values of Re x 
other than 10 8 , the same approximate procedure 
as recommended for Ree as the variable can be 
employed. For this case, however, i stea f 


figures 7(a) and (c), figure 14 and Eckert's refer- 
ence temperature should be utilized. 

To obtain an approximate relation between 
Stanton number and Re x the following procedure 
is recommended : Find the value of Re$ corre- 
sponding to the specified Re x by interpolation or 
extrapolation of figure 12. From this value of 
R(e find the Reynolds analogy factor by similar 
use of figure 11 for the specified values of f t Jt aw 
and A/$. This value of the Reynolds analogy 
factor and the value of C/ obtained as shown in 
the previous paragraph are sufficient to solve for 
the Stanton number for the specified conditions. 

CLOSING REMARKS 

No attempt has been made in this analysis to 
include the effects of dissociation, shock waves, 
radiation, slip flow, or induced pressure gradients. 

A rough estimate of the ('fleet of dissociation 
may be inferred, as pointed out by Eckert (ref. 3), 
from the theory of laminar boundary layers. 
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Tims, both Crown (ref. 20) and Moore (ref. 21) 
conclude that the effect of dissociation on friction 
factor and heat flux for the laminar boundary 
layer will be small if the wall temperature is less 
than the air dissociation temperature, which, even 
at a pressure of 0.0001 atmosphere, is above 
3000° F. Their calculations were made for Mach 
numbers up to 20. Where dissociation is appre- 
ciable, it is recommended that the heat-transfer 
coefficients presented in this report be interpreted 
as based on an enthalpy difference instead of a 
temperature difference. Thus, 

A=_*e2_ 

H»-n aw 

and 


St = 


<1 

ijj w 7/d ip) 'W sP < 5 


where II is the enthalpy. 

Although in practice there would be a shock 
wave originating near the leading edge of the 
flat plate for high Mach numbers, the effect on 
temperature and pressure distributions appears 
too complicated to be taken into account. There- 
fore, constancy of free-stream pressure and tem- 
perature has been assumed. 

The possibility of encountering slip flow at high 
Mach number must also be considered. Accord- 
ing to Eckert (ref. 3), the assumption of a con- 
tinuum is valid as long as the Knudsen number 


h 


77 r 


Ms 


-.Vs 

0.499 Ret} 


is loss than 0.01. To obtain a con- 


servative estimate for the range of conditions 
considered in this report, values of A/* and Re $ 
of 20 and 10 4 , respectively, are used. For these 


values the Knudsen number is 0.00298, which is 
well below the criterion for slip flow. 

SUMMARY OF RESULTS 

The following results were obtained from the 
analysis of turbulent flow and heat transfer over 
a flat plate at high Mach numbers: 

1. The frictional heating that occurs at high 
Mach numbers produced a flattening of the 
velocity profile, as does heating the plate by other 
means. Cooling the plate caused the velocity 
gradients near the outer edge of the boundary 
layer to increase. 

2. The skin-friction coefficients and Stanton 
numbers at a given Reynolds number decreased 
as Mach number increased. 

3. The curves for t lie ratio of friction coefficient 
to the incompressible coefficient against Mach 
number agreed closely with the curves for the 
ratio of Stanton number to incompressible Stanton 
number against Mach number. 

4. Cooling the plate to offset the effects of 
frictional heating increased the friction coefficients 
and Stanton numbers. 

5. Frictional heating at high Mach numbers 
produced a thinning of the boundary layer at a 
given position on the plate for the same free- 
stream properties. 

(k The predicted friction coefficients and Stan- 
ton numbers agreed closely with representative 
experimental data. 

7. The Reynolds number effect on both friction 
factor and Stanton number increases greatly with 
increasing Mach number. 

Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, January 17, 11)58 



APPENDIX A 

DERIVATION OF TURBULENT MOMENTUM AND ENERGY EQUATIONS 


MOMENTUM EQUATION 

The momentum equation for compressible 
boundary-layer flow past a flat plate can be written 
as 


du . bu b ( bu\ 
pU bx pl by by \ by) 

and the continuity equation as 

b(pu) dfp?Q .. 
bx ^ by 


(Al) 


(A2) 


from assuming that the laminar and turbulent 
shear stresses are of the same order of magnitude. 
The seventh is consistent with the sixth, since a 
triple correlation should be roughly of the mag- 
nitude of a double correlation raised to the 3/2 
power The eighth appears justified since it might 
be expected that k* and p' should be at least one- 
half order of magnitude less than k and p. 

With the preceding criteria, the time-averaged 
momentum equation becomes, on neglecting terms 
of magnitude S and less, 


Time derviatives and pressure gradients are 
neglected in equations (Al) and (A2), as they drop 
out when time averages are taken. 

The instantaneous quantities in equation (Al) 
are now replaced by their time averages and 
fluctuating components, which are written as 


u=u-j-u' 


-p+p'\ 

^p+p'/ 


p=p+p'’ 


(A 3) 


and time averages are taken term by term. The 
following order-of-magnitude criteria are used for 
both the momentum and energy equations: 


dr 


*0(1) 


p f u,t «0(1) 
r»0 (<$) 

0(5 2 ) 

u'r' , pV, etc. «0(5) 


P ' u' v ' , p ■' u ' 2 , e tc . ~ 0 (5 3/2 ) 

Double correlations containing k' and p'«0(6 :i ) 

The first five of these criteria are the usual 
boundary-layer assumptions. The sixth results 
28 


— bu — ~~TT\ du b ( - bu — ; — y 

pu - +CPV+P V ) ^— dy (p ^-pU V , 

and tl e continuity equation, 

bCpu) . d(p?) , 5(pV) 


by 


=0 


bx by 

Considering the relation 

pF+p V = pr 

equations (A4) and (A5) can be rewritten 
— bu . — bu b / _ du --r-r N 

nit — UrtJ' — I ii nit'iV 


(A4) 


(A5) 


(A6) 


~ pu 6J +pr dy^dy pU '° ') (A7) 


and 


dipu) . d(pv) 


dx 


dy 


(AS) 


Comparison of equations (Al) and (A7) leads to 
the definition of r as 


_ du — — 

T = jU — — p u v 

dy 


ENERGY EQUATION 


(A9) 


Tin* energy equation for compressible boundary- 
layer Jlow past a flat plate is 


dt , 

C'C" 0j . + P' r r . 


at_ d / dt\ /d«\* A1 _, 

^dJj-dyVdyJ+^Kdy) (Al0) 
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Time derivatives and pressure gradients are 
again neglected, as time-averaging cancels them. 
If the momentum equation (Al) is multiplied 
through by u and then added to equation (A10), 
the result is 


P u 




(All) 


where c v is considered constant. 

Again substituting for the instantaneous quan- 
tities the sum of the time-averaged ami fluctuating 
components, and then neglecting terms of small 
order of magnitude on taking time averages, 
equation (All) becomes 

~ pTl ir ( C ” 1+ Y) + ( ^ +57F) iy (^i) 

L % (* ^-- P c p —-pu—) (A12) 


Again employing equation (A6), equation (Al2) 
becomes 

-pu ^ (V+y)+^ (m+y) 

(k + ~—pc p v 't' — p uu> v f ^ (A 13) 

dy\ dy dy / 


A comparison of equations (All) and (Al 3) leads 
to the definition of q as 


_ du 

u 

dy 

It should be noted that this treatment gives 
no density fluctuation terms in the expressions 
for r and q. This same result was found by 
Van Driest (ref. 22) and by Rubesin (ref. 19). 
Combining pV with pv and writing the sum as 
pv present no difficulty, because, in a complete 
solution, pv could be eliminated from the momen- 
tum and energy equations by t he e quation of 
continuity. An assumption for pV would be 
necessary only if it were desired to calculate r. 



pCpV't'— puu'v' ) (A14) 


) 



APPENDIX B 


PROOF THAT — = FOR Pr-a = 1 


In terms of eddy diffusivities, the momentum 
and energy equations may be written as follows: 


— du , — du d r , _ 

pw d7 +pi, ^ = d^L (#1+pe 


■> I ] (B1) 


( rJ+ 'i) +p,, §y( C » J+ Y) 


b 

by 


\ [( E +^<») ^+m(m+P*) §£] (B2) 

The energy equation (B2) can be rearranged to 
read 


h ( r " 7+ i>) +p '’ § y ( r ” ,+! f) 


~ ^.+ p <A b7, ) + < M + P«) ^ (y) (BH) 


-tyL7v™a* kvj-hp+p.) ^ 

For Pr=a= 1 , equation (B3) becomes 
p " dr ( r " 7+ 2' ) +p? b7, ( ( '" 7 + t) 

= |/ [ (p+p€) dT/ ( r,,7+ y)] (B4) 

If equations (Bl) and (B4) are each solved for 
pt 7 and the results equated, there is obtained 


The expression given by equation (B6) for 

— 

Cpt + 7 j- is substituted into the right side of equa- 
tion (Bo). Since the right side becomes identical 
with the left side upon this substitution, equation 
(B6) is a valid relation. 

The constants A and B in equation (B6) are 
evaluated as follows: At y= 0, 

w 0, t , B — Cpt if 


-du_ T dt q w 

P by Tu ” k by iw ’ A ~ r. 

With the constants thus evaluated, equation (B6) 
becomes (dropping the bars for convenience) 


c,t+\=~-u-Vc v t v (B7) 

" T w 


If (‘qua ion (B7) is made dimensionless, it becomes 
simply 

u + = 7 M (B8) 


Previoi sly obtained relations for t/t 1p and 


are 


a r 

,-i-i b»/ 

— bu 

T 

by L 

iM-f-pe) 

by_ 

pl ' by 

T w 


du 

6jc 


d 


_ ! p - p<i ct(^ 7+ f )]~ pp cK ^+f) 


J ■_( p t_P_ ( \ * ^ + 

r i« VPic Pw ptclptc/ dy+ 

( k P t \ dt + 

\k, Pr„p w a ix w / Pw ) dy + 


(JiQw 


(5) 




The assumption is now made that 


e,l+%=Au+H 


(B5) 

(B6) 


, « + / M , P « \ d^ + , n 

0 Pw Pw/PwJ d?/ + ^ 


If use is made of equations (B8), (10), and (26) 
and the fact that Pr=a= 1, the resulting equation 
can be educed to 


(B9) 


30 



APPENDIX C 


VELOCITY AND TEMPERATURE PROFILES FOR CONSTANT RATIO OF TURBULENT TO MOLECULAR 
SHEAR STRESS AT yf AND WITH MOLECULAR SHEAR STRESS AND HEAT TRANSFER INCLUDED 

IN REGION AWAY FROM WALL (Pr = l) 


The equation for velocity profile; used near the 
wall is equation (16), where 1 — at + ' = 1 — au +2 for 
Pr=l and 0= 0 (see eqs. (BS), (26), and defini- 
tions of P and t + '). 

Previously, the expression for r/r M was shown 
to be 

r u d u + . p € d u + /-x 

— -T- + + - T+ W 

t w ay* p w Pic/Pw u y 

and far from the wall (eq. (8)), 



Vdy T V 


Combining equations (5) and (8) and assuming 
constant shear stress across the boundary layer 
give 



Substitution of equations (02) and (03) into (Cl) 
yields 


1 = 


/ d?/ + 


(1 -aw +2 ) 0 ' 6s + 


l 

1 —au +2 / d 2 u + \ 2 I 

\d // + 2 / J 


k 2 Vd ?/ 


d u + 

d y + 


(C4) 


Solving for d 2 u + /dy+ 2 gives 

du + ^ 2 


d 2 u + 

dy +2 




)• 


■^(l-ortC*) [l-(l -«»+’)»•<* 


(Co) 


If a change in variables is made as 

d-M + 

equation (05) can be integrated to give 




J >W+ K(]U + 

Y ( 1 — a U *'2 ) [ 1 — ( 1 


— a— (0(3) 


The variations 
temperature are 


of density and 





P Or 

Pto t- 


viscosity with 


(C2) 


For Pr—a— 1 
expressed as 


and 0=0, the temperature ratio is 


1 = l-«»+* (C3) 

* tc 


The solutions for as a function of y + can be 
obtained by a process of iteration. Assumed 
values of v for a given increment in u + are sub- 
stituted into the right side of equation (C6) until 
it equals the left side. The relation between w + 
and ?yO is then calculated from 





d?/ + 

v 


(C7) 


From equation (B8), u + = T+ for qlq v =TjT Vf so that 
the relation between and y + is also known. 
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APPENDIX D 


VELOCITY PROFILES FOR LINEAR VARIATION OF SHEAR STRESS AND ENERGY TRANSFER 

ACROSS BOUNDARY LAYER (Pr=l) 


Because the variation of shear in the thin region 
near the wall is negligible, the same equations are 
used in the present ease as were used in appendix 
C. 

From equation (5), the equation for r/r w far 
from the wall, neglecting the viscous stress, is 


r p e dw + 

r to P to M to/p to d ^ 

The expression for e/(pjp w ) far from the wall is, 
from equation (8), 



For a linear variation in shear stress, 


-= 1 - 


y 

5+ 


(Dl) 


Combining the foregoing equations yields 



( 1 ) 2 ) 


For 0= 0, Pr=a= 1 (r/r w =^ w ), 


/-=l-orte+*=^ 

*w p 


(D3) 


Substituting equation (D3) into (D2), rearrang- 
ing, and taking the square root of both sides give 


d 2 ?/ 1 
dy +2 


/ dtt + V 

— * vh/ + / 


V 1 — OLU^ 


V-4: 


+ 

<5~ + 


(D4) 


Letting v—du' i '/dy + in equation (1)4) and inte- 
grating give 


r=Pj e 


k Cw+ d u+ 

J “t y'i-Qid 2 /i--4_ 

A 


-± f u 

* + J«t 


tt+ du+ 


(1)5) 


Equation (Do) can be solved by iteration for 0=0 
and a given a and <5 + , to give ?/+ as a function of 
?/ + - From equation (118), u + =T + for q!q w — t/t*, 
so that the relation between Y T+ and ?/*' is also 
known. 


REFERENCES 


1* Rubesin, Morris \\ .. Maydew, Randall C\, and Varga, 
St(‘vt‘n A.: An Analytical and Experimental In- 
vestigation of the Skin Friction of the Turbulent. 
Boundary Layer on a Flat Plate at Supersonic 
Speeds. X A C A T X 230 5 , 105!. 

2. ( oles, Donald: Measurements in the Boundary Layer 

on a Smooth Flat Philo in Supersonic Flow. Ill — . 
Measurements in a Flat -Plate Boundary Layer at 
the Jet Propulsion Laboratory. Rep. Xo. 20 71, 
Jet Prop. Lab., (M.T., June 1, 1053. (Contract 
Xo. DA 01 105 OKI) LS, Dept. Army, Ord. Corps.) 

3. Eckert, Ernst R. (L: Survey on Heat Transfer at 

High Speeds. Tech. Rep. 54 70, Aero. Res. Lab., 
\\AI)(\ Apr. 1054. (Contract Xo. AF 33(616)- 
22 4 .) 

4. Deissler, Robert ( 5 . : Analytical and Experimental In- 

vestigation of Adiabatic Turbulent Flow in Smooth 
Tubes. XAC A TX 2138, 1050. 

5. Deissler, R. (L: Invest igat ion of Turbulent Flow and 

32 


He; t Transfer in Smooth Tubes, Including the 
Elf ‘Clsof \ ariable Fluid Properties. Trans. ASM K, 
vol 73, no. 2, Feb. 1051, pp. J01 107. 

0. Deiss er, R. (r.: Heat Transfer and Fluid Friction for 
Ful y Developed Turbulent. Flow of Air and Super- 
crit.eal Water with Variable Fluid Properties 
Trans. ASM E, vol. 7(>, no. 1, Jan. 1054, pp. 73 So’ 

7. von Kilrman, Th.: Turbulence and Skin Friction. 

Jou ■. Aen>. Sci., vol. J, no. 1, Jan. 1034, pj). 1 20. 

8. Deissler, R. fb, and Eian, C. S.: Analytical and Ex- 

perimental Investigation of Fully Developed Tur- 
buh nt Flow of Air in a Smooth Tula 1 with Heat 
Transfer with Variable Fluid Properties. XAC A 
TX 2020, 1052. 

0. Deissl *r, Robert ( L: Turbulent Heat. Transfer and 
l 1 rid ion in the Entrance Regions of Smoot h Pas- 
sages. Trans. ASME, vol. 77, no. 8, Xov. 1955, 
PP* 1221-1232; discussion, p. 1233. 




ANALYSIS OF TURBULENT FLOW AND HEAT TRANSFER ON FLAT PLATE AT HIGH MACH NUMBERS 33 


10. Deissler, Robert G.: Analysis of Turbulent Heat 

Transfer, Mass Transfer, and Friction in Smooth 
Tubes at High Prandtl and Schmidt Numbers. 
NACA Hep. 1210, 1955. (Supersedes NAG A TN 
3145.) 

11. Deissler, R. G., and Loeffler, A. L., Jr.: Turbulent 

Flow and Heat Transfer on a Flat Plate at High 
Mach Number with Variable Fluid Properties. 
Paper No. 55-A-133, AS ME, 1955. 

12. Reichardt, H.: Fundamentals of Turbulent Heat 

Transfer. Arehiv f. d. gesamte \\ armetech., no. 
6/7, 1951, pp. 129 -142. 

13. Goldstein, S., ed.: Modern Development in Fluid 

Dynamics. Vol. II. Clarendon Press (Oxford), 
1938, p. 351. 

14. Laufer, John: The Structure of Turbulence in Fully 

Developed Pipe Flow. NACA Rep. 1174, 1954. 
(Supersedes NACA TN 2954.) 

15. KlebanofT, P. S.: Characteristics of Turbulence in a 

Boundary Layer with Zero Pressure Gradient. 
NACA Rep. 1247, 1955. (Supersedes NACA TN 
3178.) 

16. von Karm&n, Th.: The Problem of Resistance in Com- 

pressible Fluids. Reale Accad. d’ltaliu (Rome), T. 
XIII, Sept,- Oct.. 1935. 

17. Monaghan, R. J., and Cooke, J. R.: 1 he Measurement 

of Ileat Transfer and Skin Friction at Supersonic 
Speeds. Pt. IV. Tests on a Flat Plate at M =2.82. 
Tech. Note AFRO. 2171, British RAF, June 1952. 

18. Hill, F. K.: Boundary-Layer Measurements in Hyper- 

sonic Flow. Jour. Aero. Sci., vol. 23, no. 3, Jan. 
1956, pp. 35-42. 

19. Rubesin, Morris W.: A Modified Reynolds Analogy 

for the Compressible Turbulent Boundary Layer on 
a Flat Plate. NACA TN 2917, 1953. 

20. Crown, J. C.: The Laminar Boundary Layer at Hy- 

personic Speeds. N A VOH I) Rep. 2299, Naval Ord. 
Lab., Apr. 15, 1952. 

21. Moore, L. L.: A Solution of the Laminar Boundary- 

Layer Equations for a Compressible Fluid with 
Variable Properties, Including Dissociation. Jour. 
Aero. Sci., vol. 19, no. 8, Aug. 1952, pp. 505 5 IS. 

22. Van Driest, E. R.: Turbulent. Boundary Layer in Com- 

pressible Fluids. Jour. Aero. Sci., vol. 18, no. 3, 
Mar. 1951, pp. 145 161. 


23. Lobb, It. Kenneth, Winkler, Eva M., and Persh, 

Jerome: Experimental Investigation of Turbulent 
Boundary Layers in Hypersonic Flow. Jour. Aero. 
Sci., vol. 22, no. 1, Jan. 1955, pp. 1-9; 50. 

24. Brinich, Paul F., and Diaconis, Nick S.: Boundary- 

Layer Development and Skin Friction at Mach 
Number 3.05. NACA TN 2742, 1952. 

25. Spivaek, H. M.: Experiments in the Turbulent 

Boundary Layer of a Supersonic Flow. Rep. CM- 
615 (AL 1052), Aerophys. Lab., North Am. Avia- 
tion, Inc., Jan. 16, 1950. 

26. Monaghan, R. J., and Johnson, J. E.: The Measure- 

ment of Heat Transfer and Skin Friction at Super- 
sonic Speeds. Pt. II. Boundary-Layer Measure- 
ments on a Flat Plate at Af= 2.5 and Zero Heat 
Transfer. C.P. 64, British ARC, 1952. 

27. Schultz-Grunow, F.: New Frictional Resistance Law 

for Smooth Plates. NACA TM 986, 1941. 

28. Slack, Fllis G.: Experimental Investigation of Heat 

Transfer Through Laminar and Turbulent Boundary 
Layers on a Cooled Flat Plate at a Mach Number 
of 2.4. NACA TN 2686, 1952. 

29. Pappas, C. Measurement of Heat Transfer in the 

Turbulent Boundary Layer on a Flat Plate in 
Supersonic Flow and Comparison with Skin-Friction 
Results. NACA TN 3222, 1954. 

30. Emmons, M. A., Jr., and Blanchard, R. F.: Prelimi- 

nary Investigation of the Transfer of Heat from a 
Flat Plate at a Mach Number of 1.5. NACA RM 
L511131, 1951. 

31. Jacob, Max, and Dow, W. M.: Heat. Transfer from a 

Cylindrical Surface to Air in Parallel Flow with 
and without Unheated Starting Sections. Trans. 
ASM F, vol. 68, no. 2, Feb. 1946, pp. 123 134. 

32. Dhawan, Satish: Direct Measurements of Skin Fric- 

tion. NACA Rep. 1121, 1953. (Supersedes NACA 
TN 2567.) 

33. Kempf, Gunther: Neue Frgebnisse der \\ iderstands- 

forschung. Werft Reederei Hafen, Bd. 10, June 
1929, pp. 234 239; 247 253. 

34. Trass, Olev: Sublimation Mass Transfer Through 

Compressible Boundary Layers on a Flat Plate. 
Ph. D. Thesis, M.T.T., Sept. 1958. 


U.S. GOVERNMENT PRINTING OFFICE: I960 




